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£<\ Abstract 

We consider the effect of the wind and the dissipation on the nonlinear 
stages of the modulational instability. By applying a suitable transforma- 

^"-"l tion, we map the forced/damped Nonlinear Schrodinger (NLS) equation into 

the standard NLS with constant coefficients. The transformation is valid 
as long as \Tt\ <C 1, with V the growth/damping rate of the waves due to 
the wind/dissipation. Approximate rogue wave solutions of the equation are 

^h presented and discussed. The results shed some lights on the effects of wind 

and dissipation on the formation of rogue waves. 
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1. Introduction 



Modulational instability, also known as the Benjamin- Feir instability in 
J> the water wave community, has been discovered in the late sixties indepen- 

dently by Benjamin and Feir [TJ and Zakharov [2] (see [3] for an historical 
review on the subject and possible applications). It describes the exponen- 
tial growth of an initially sinusoidal long wave perturbation of a plane wave 
solution of the one dimensional water wave problem. For water waves the 
condition of instability in infinite water depth is that 2\/2a k N > 1, where 
a Q is the amplitude of the plane wave and k is the corresponding wave 
number; N = ko/AK is the number of waves under the perturbation of 
wavenumber AK. The modulaitonal instability is frequently studied within 
the Nonlinear Schrodinger (NLS) equation that describes weakly nonlinear 
and dispersive waves in the narrow band approximation. In this context, 
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the nonlinear stages of the modulational instability are described by exact 
solutions of the NLS, known as Akhmediev breathers [U |5]. Other exact 
NLS solutions which describe the focussing of an initially non-small pertur- 
bation have been derived in [61 [7] . Such solutions have been considered as 
prototypes of rogue waves [HI E] ■ 

Within the one dimensional NLS equation, the modulational instability is 
well understood. What is probably less clear is the modulation of waves and 
the formation of rogue waves in forced (by wind) or damped (by dissipation) 
conditions. In this regard in the past there has been a number of experimental 
works, PUJ El H2] , which did not gave a clear picture on the effect of the wind 
on the modulaitonal instability. A careful discussion of the discrepancy of the 
results presented in the above papers can be found in [12J. According to their 
discussion the role of the wind is twofold: i) the wind changes the growth 
rate of the instability; ii) the natural selection of the sideband frequency is 
altered with respect to the no wind conditions. 

Concerning damping effects, it has been showed in [T3] that any amount 
of dissipation stabilizes the modulational instability, questioning the role of 
the modulational instability in the formation of rogue waves, |14j . More 
recently, the role of dissipation and wind in the modulational instability has 
been considered together within the NLS equation, [15] (then confirmed by 
fully nonlinear simulations, [IS])- The authors performed a linear stability 
analysis and numerical simulations and found that, in the presence of wind, 
young waves are more sensitive to modulational instability than old waves. 

The just mentioned numerical results (except the one in [16]) are all based 
on the following forced and damped Nonlinear Schrodinger equation: 

'£-4?-^=™- (i) 

A is the wave envelope, a and /3 are two coefficients that depend on the 
wavenumber, ko, of the carrier wave. The right-hand side is responsible for 
the forcing, T > 0, and/or dissipation, T < 0. The two effects are additive 
so that T is in general the sum of forcing coefficient plus a damping one. 
The wind forcing depends on the ratio between air and water density and 
the dissipation on the water viscosity, therefore the absolute value of T is 
always a small quantity. Finding analytical solutions of equation (flj is not 
an obvious task. In the present paper we take advantage of the smallness of T 
and, after a suitable transformation, we are able to find breather solutions of 



the forced-damped NLS equation. In the following sections we first describe 
the transformation and then present the rogue wave analytical solutions. 

2. Reduction of the forced/damped NLS to the standard NLS 

We considered the NLS equation discussed in [T5] 



with 
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here k is the Von Karman constant and u* is the friction velocity, g is the 
gravity acceleration, p a and p w are the air and water density, respectively; 
7 is a coefficient to be determined from the solution of Rayleigh equation 
associated to the stability of the wind wave problem (see also [17] for a 
justification of the wind forcing term); c is the phase velocity, v is the water 
kinematic viscosity. In [T5] the equation is written in a nondimensional form 
and the coefficient K = T/u is introduced). The surface elevation is related 
to the envelope as follows: 

r)(x, t) = -i A(x, t) exp[i(k x - u t)} + c.cj . (4) 

Note that we use a different definition of the surface elevation from the one 
in [15] where the 1/2 factor is not included (the consequence is that the 
coefficient in the nonlinear term in equation (J2J) differs by a factor of 4 from 
the one in equation (3.1) in [15]). If e is the small parameter in the derivation 
of the NLS, then it is assumed that the right-hand side term in pi) is of the 
order of e 2 as the nonlinear and the dispersive term. 
We consider the following new variable: 

B(x,t) = A(x,t)e- rt (5) 

and by selecting a coordinate system moving with the group velocity we get: 

^f -a|f-/3exp-W5 = (6) 



were a and are the coefficients of the dispersive and nonlinear term, re- 
spectively. Written in the above form the effect of the forcing/damping term 
enters as a factor in front of nonlinear term and has the role of enhanc- 
ing/decreasing the nonlinearity of the system as the wave evolve in time. 
Recalling that T is usually small, we Taylor expand the exponential and 
re-write the equation as follows: 

i^-ay-Mt)\B\ 2 B^0 (7) 

with p(t) = 1/(1 — 2Tt). Let's introduce the following change of coordinates: 

X {x,t)=p(t)x, T{t)=p(t)t (8) 

and scale the wave envelope function B as follows 
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After the transformation, the equation rt6J) results in: 



^ ^ ^-^-m 2 ^P = (10) 

i.e., the NLS equation with constant coefficients. We have transformed the 
forced/damped NLS equation into the standard NLS equation whose solu- 
tions can be studied analytically. From a physical point of view the trans- 
formation (and consequently the validity of the solutions) is valid as long 
as 2|Tt| <C 1 (the transformation is singular for 2|Tt| = 1). We underline 
that the transformation of the forced/damped NLS equation to the standard 
one has been possible only for l/p(t) equal to a linear function in t. For 
other functional dependences, the transformation does not seem to be possi- 
ble. Our result is consistent with ones reported in [HI [19] where analytical 
solutions of the variable coefficient NLS equation are described. 

3. Rogue wave solutions 

In the following we will present three analytical solutions corresponding 
to the Peregrine, the Akhmediev and the Kuznetsov-Ma breathers for the 
standard NLS. 



The Peregrine solution also known as rational solution, has been origi- 
nally proposed in [5] . It has the peculiarity of being not periodic in time and 
in space: it is a wave that "appears out of nowhere and disappears without 
trace" [20l [2Tj ; its maximum amplitude reaches three times the amplitude of 
the unperturbed waves. For the above reasons it has been considered as spe- 
cial prototype of freak wave, [2T]. The Peregrine solution has been recently 
reproduced experimentally in wave tank laboratories [22] and in optical fibers 
|23j . Below we present an exact analytical solution of equation ([7]) which is 
the analogous of the Peregrine solutions but for the forced/damped case: 

B(xt)-BC(xt)( 4(1 - »2/?Bgp(t)t) \ 

B(x, t) - B G(x, t) { a + a{2/3B 2 p{tW + 2l3B 2 Mt)x) 2 ~ l ) C 11 ) 



(12) 



with 

G(x,t) = v W)ew i (p^- - PBlp{t)t\ 

In figure [T] we show an example of such solution for steepness 0.1 and forcing 
coefficient K = 0.0004 (the same value has been used in [15]). The axis are 
normalized by the wave period, the wavelength and the initial wave amplitude 
B . The effect of the wind/dissipation is to increase/reduced the amplitude 
of the plane wave. As in the case of the standard NLS, the wave appears 
only once in time and space. 

The Akhmediev solution jlj describes the modulational instability in its 
nonlinear regime; it is periodic in space. It is characterized by an amplifi- 
cation factor which ranges from 1 to 3 (this last value corresponds to the 
Peregrine solution). In the presence of a forcing or damping, the breather 
has the following analytical form: 

B (x, t) = b oG (x, t) ( )^M^m^^^^jm\ _ \ (13) 

K J \V2cosh[ap{t)t]-V2-u 2 cos{up{t)x} j 

and 

u=%, * = ^/f, ~° = vy[2^V\ a = (3B 2 a. (14) 

The function G(x,t) is reported (12). It should be noted that the function 



is periodic in space with a period that changes in time. In figure [2] we show 
an example of such solution for steepness 0.1, iV = 5 and forcing coefficient 
oiK = 0.0004. 
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Figure 1: The Peregrine solution of the forced NLS equation. 
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Figure 2: The Akhmediev solution of the forced NLS equation. 
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Figure 3: The Kuznetsov-Ma solution of the forced NLS equation. 



The Kuznetsov-Ma solution [H] is periodic in time and decrease exponen- 
tially in space. While for the Akhmediev breather the large time (positive or 
negative) limit is a plane wave plus a small perturbation, the modulation for 
the Ma breather is never small. The solution for the forced/damped equation 
is here reported: 



B(x,t) = B G(x,t) 



-v^/i 2 cos[pp(t)t] + iy/2p sin[pp(t)t)} 
y/2 cos[pp(t)t] — a/2 + jl 2 cosh.[vp(t)x] 



with 



/' 




P 



pV2TW, P = PB*p. 



(15) 



(16) 



ft is a parameter related to the amplification factor. In figure [3] we show an 
example of such solution for steepness 0.1, jl = y/2 and forcing coefficient of 
K = 0.0004. The periodicity (appearance of maxima) changes in time and 
increase in the presence of forcing and decrease for the damping case. 



4. Conclusion 

In the present Letter we have considered the problem of generation of 
rogue waves in the presence of wind forcing or dissipation. Our work is 
based on the one dimensional forced/damped NLS equation. 



Under the assumption of 2\Tt\ <C 1, where T is the forcing (T > 0) or 
damping (T < 0) term, we have shown how the equation can be mapped in 
the standard NLS equation with constant coefficients, fn this framework, we 
have found explicit analytical breather solutions. 

As mentioned the effect of wind/dissipation is to increase/reduce in time 
the coefficient in front of the nonlinear term. This has an impact on the 
modulational instability; in particular, an initially stable (unstable) wave 
packet could be destabilized (stabilized) by the wind (dissipation). Similar 
results have been obtained for the interaction of waves and current (see [24"]). 
The present results should be tested in wind waves tank facilities. 
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